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ABSTRACT 

We present the application of the Fast Independent Component Analysis (FastICA) 
technique for blind component separation to polarized astrophysical emission. We 
study how the Cosmic Microwave Background (CMB) polarized signal, consisting 
of E and B modes, can be extracted from maps affected by substantial contamina- 
tion from diffuse Galactic foreground emission and instrumental noise. We implement 
Monte Carlo chains varying the CMB and noise realizations in order to asses the av- 
erage capabilities of the algorithm and their variance. We perform the analysis of all 
sky maps simulated according to the Planck satellite capabilities, modelling the sky 
signal as a superposition of the CMB and of the existing simulated polarization tem- 
plates of Galactic synchrotron. Our results indicate that the angular power spectrum 
of CMB £'-mode can be recovered on all scales up to £ ~ 1000, corresponding to the 
fourth acoustic oscillation, while the i3-mode power spectrum can be detected, up to 
its turnover at £ ~ 100, if the ratio of tensor to scalar contributions to the temperature 
quadrupole exceeds 30%. The power spectrum of the cross correlation between total 
intensity and polarization, TE, can be recovered up to £ ~ 1200, corresponding to the 
seventh TE acoustic oscillation. 

Key words: methods - data analysis - techniques: image processing - cosmic mi- 
crowave background. 



1 INTRODUCTION 

We are right now in the epoch in which the cosmological observations are revealing the finest structures in the Cosmic 
Microwave Background (CMB) anisotropies^ . After the first discovery of CMB total intensity fiuctuations as measured by the 
COsmic Background Explorer (COBE) satellite (see Smoot 1999 and references therein), several balloon-borne and ground- 
based operating experiments were successful in detecting CMB anisotropies on degree and sub-degree angular scales (De 
Bernardis et al 2002, Halverson et al. 2002, Lee et al. 2001, Padin et al. 2001, see also Hu & Dodelson 2002 and references 
therein). The Wilkinson Microwave Anisotropy Probe (WMAP, see Bennett et al. 2003a) satellite^ released the first year, 
all sky CMB observations mapping anisotropies down to an angular scale of about 16' in total intensity and its correlation 
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with polarization, on five frequency cliannels extending from 22 to 90 GHz. In the future, balloon-borne and ground based 
observations will attempt to measure the CMB polarization on sky patches (see Kovac et al. 2002 for a first detection); the 
Planck^ satellite, scheduled for launch in 2007 (Mandolcsi ot al. 1998, Puget et al. 1998), will provide total intensity and 
polarization full sky maps of CMB anisotropy with resolution <^ 5' and a sensitivity of a few jjK, on nine frequencies in the 

range 30-857 GHz. A future satellite mission for polarization is currently under study''. 

Correspondingly, the data analysis science faces entirely new and challenging issues in order to handle the amount of 
incoming data, with the aim of extracting all the relevant physical information about the cosmological signal and the other 
astrophysical emissions, coming from extra-galactic sources as well as from our own Galaxy. The sum of these foreground 
emissions, in total intensity, is minimum at about 70 GHz, according to the first year WMAP data (Bennett et al. 2003b). In 
the following we refer to low and high frequencies meaning the ranges below and above that of minimum foreground emission. 

At low frequencies the main Galactic foregrounds arc synchrotron (see Haslam et al. 1982 for an all sky template at 408 
MHz) and free-free (traced by Ha emission, sec Haffner, Reynolds, & Tufte 1999, Finkbeiner 2003 and references therein) 
emissions, as confirmed by the WMAP observations (Bennett et al. 2003b). At high frequencies. Galactic emission is expected 
to be dominated by thermal dust (Schlegel, Finkbeiner, Davies 1998, Finkbeiner, Schlegel, Davies 1999). Moreover, several 
populations of extra-galactic sources, with different spectral behavior, show up at all the frequencies, including radio sources 
and dusty galaxies (see Toffolatti et al. 1998), and the Sunyaev-Zeldovich effect from clusters of galax;ies (Moscardini et al. 
2002). Since the various emission mechanisms have generally different frequency dependencies, it is conceivable to combine 
multi-frequency maps in order to separate them. 

A lot of work has been recently dedicated to provide algorithms devoted to the component separation task, exploiting 
different ideas and tools from signal processing science. Such algorithms generally deal separately with point-like objects like 
extra-galactic sources (Tenorio et al. 1999, Vielva et al. 2001), and diffuse emissions from our own Galaxy. In this work we 
focus on techniques developed to handle diffuse emissions; such techniques can be broadly classified in two main categories. 

The "non-blind" approach consists in assuming priors on the signals to recover, on their spatial pattern and frequency 
scalings, in order to regularize the inverse filtering going from the noisy, multi-frequency data to the separated components. 
Wiener filtering (WF, Tegmark, Efstathiou 1996, Bouchet, Prunet, Sethi 1999) and Maximum Entropy Method (MEM, Hobson 
et al. 1998) have been tested with good results, even if applied to the whole sky (Stolyarov et al. 2002). Part of the priors can 
be obtained from complementary observations, and the remaining ones have to be guessed. The WMAP group (Bennett et 
al. 2003b) exploited the available templates mentioned above as priors for a successful MEM-based component separation. 

The "blind" approach consists instead in performing separation by only assuming the statistical independence of the 
signals to recover, without priors cither for their frequency scalings, or for their spatial statistics. This is possible by means of 
a novel technique in signal processing science, the Independent Component Analysis (ICA, see Amari & Chichocki 1998 and 
references therein). The first astrophysical application of this technique (Baccigalupi et al. 2000) exploited an adaptive (i.e. 
capable of self-adjusting on time streams with varying signals) ICA algorithm, working successfully on limited sky patches for 
ideal noiseless data. Maino et al. (2002) implemented a fast, non-adaptive version of such algorithm (FastICA, see Hyvarinen 
1999) which was successful in reaching separation of CMB and foregrounds for several combinations of simulated all sky 
maps in conditions corresponding to the nominal performances of Planck, for total intensity measurements. Recently, Maino 
et al. (2003) were able to reproduce the main scientific results out of the COBE data exploiting the FastICA technique. 
The blind techniques for component separation represent the most unbiased approach, since they only assume the statistical 
independence between cosmological and foreground emissions. Thus they not only provide an independent check on the results 
of non-blind separation procedures, but are likely to be the only viable way to go when the foreground contamination is poorly 
known. 

In this work we apply the FastICA technique to astrophysical polarized emission. CMB polarization is expected to 
arise from Thomson scattering of photons and electrons at decoupling. Due to the tensor nature of polarization, physical 
information is coded in a entirely different way with respect to total intensity. Cosmological perturbations may be divided 
into scalars, like density perturbations, vectors, for example vorticity, and tensors, i.e. gravitational waves (sec Kodama & 
Sasaki 1984). Total intensity CMB anisotropics simply sum up contributions from all kinds of cosmological perturbations. 
For polarization, two non-local combinations of the Stokes parameters Q and U can be built, commonly known as E and B 
modes (sec Zaldarriaga, & Seljak 1997, and Kamionkowski, Kosowsky & Stebbins 1997 featuring a different notation, namely 
gradient G for E and curl C for B). It can be shown that the E component sums up the contributions from all the three 
kinds of cosmological perturbations mentioned above, while the B modes are excited via vectors and tensors only. Also, scalar 
modes of total intensity, which we label with T in the following, are expected to be strongly correlated with E modes: indeed, 
the latter are merely excited by the quadrupole of density perturbations, coded in the total intensity of CMB photons, as seen 
from the rest frame of charged particles at last scattering (see Hu et al. 1999 and references therein). Therefore, for CMB, 

^ astro.estec.esa.nl/SA-general/Projects/Planck 
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the correlation TE between T and E modes is expected to be the strongest signal from polarization. The latter expectation 
has been confirmed by WMAP (Kogut et al. 2003) with a spectacular detection on degree and super-degree angular scales; 
moreover, a first detection of CMB E modes has been obtained (Kovac ot al. 2002). 

This phenomenology is clearly much richer with respect to total intensity, and motivated a great interest toward CMB 
polarization, not only as a new data set in addition to total intensity, but as the best potential carrier of cosmological 
information via electromagnetic waves. Unfortunately, as wo describe in the next Section, foregrounds are even less known in 
polarization than in total intensity, see De Zotti (2002) and references therein for reviews. For this reason, it is likely that 
a blind technique will be required to clean CMB polarization from contaminating foregrounds. The first goal of this work is 
to present a first implementation of the ICA techniques on polarized astrophysical maps. Second, we want to estimate the 
precision with which CMB polarized emission will be measured in the near future. We exploit the FastICA technique on low 
frequencies where some foreground model have been carried out (Giardino et al. 2002, Baccigalupi et al. 2001). 

The paper is organized as follows. In Section 2 we describe how the simulation of the synchrotron emission were obtained. 
In Section 3 we describe our approach to component separation for polarized radiation. In Section 4 we study the FastICA 
performance on our simulated sky maps. In Section 5 we apply our technique to the Planck simulated data, studying its 
capabilities for polarization measurements in presence of foreground emission. Finally, Section 6 contains the concluding 
remarks. 



2 SIMULATED POLARIZATION MAPS AT MICROWAVE FREQUENCIES 

We adopt a background cosmology close to the model which best fits the WMAP data (Spergel et al. 2003). We assume a flat 
Friedman Robertson Walker (FRW) metric with an Hubble constant Ho = lOO/i km/sec/Mpc with h = 0.7. The Cosmological 
Constant represents 70% of the critical density today, = 0.7, while the energy density in baryons is given by fli,h^ = 0.022; 
the remaining fraction is in Cold Dark Matter (CDM); we allow for a re- ionization with optical depth r = 0.05 (Becker 
et al. 2001). Note that this is a factor 2-3 lower than found in the first year WMAP data (see Bennett et al. 2003a and 
references therein), since we built our reference CMB template before the release of WMAP data. Cosmological perturbations 
are Gaussian, with spectral index for the scalar component leading to a not perfectly scale invariant spectrum, ns = 0.96, and 
including tensor perturbations giving rise to a B mode in the CMB power spectrum. We assume a ratio R = 30% between 
tensor and scalar amplitudes, and the tensor spectral index is taken to be nr = —R/6.8 according to the simplest inflationary 
models of the very early Universe (see Liddlo & Lyth 2000 and references therein). The cosmological parameters leading to 
our CMB template can be summarized as follows: 

ft = 0.7 . Oa = 0.7 , nth^ = 0.022 , Qcdm ^ 1 - Qa - ilt , r = 0.05 , ns = 0.96 , i? = 0.3 , nr = -R/6.8 . (1) 

We simulate whole sky maps of Q and U out of the theoretical and Cf coefficients as generated by CMBfast (Seljak & 
Zaldarriaga 1996), in the HEALPix environment (Gorski et al. 1999). The maps are in antenna temperature, which is obtained 
at any frequency ly multiplying the thormodynamical fluctuations by a factor cxpx/{expx— 1)'^, whore x — hu/kTcMB , h, k 
are the Planck and the Boltzmann constant, respectively, while Tcmb = 2.726 K is the CMB thormodynamical temperature. 

The polarized emission from diffuse Galactic foregrounds in the frequency range which will be covered by the Planck 
satellite is very poorly known. On the high frequency side the Galactic contribution to the polarized signal should be dominated 
by dust emission (Lazarian & Prunet 2002). The first detection of the diffuse polarized dust emission has been carried out 
recently (Benoit et al. 2003), and indicates a 3 — 5% polarization on large angular scales and at low Galactic latitudes. On 
the low frequency side, the dominant diffuse polarized emission is Galactic synchrotron. Observations in the radio band cover 
about half of the sky at degree resolution (Brouw & Spoelstra 1976), and limited regions at low and medium Galactic latitudes 
with 10' resolution (Duncan et al. 1997, Uyaniker et al. 1999, Duncan et al. 1999). Analyses of the angular power spectrum of 
polarized synchrotron omission have boon carried out by several authors (Tucci et al. 2000, Baccigalupi et al. 2001, Giardino 
et al. 2002, Tucci et al. 2002, Bernardi et al. 2003). 

Polarized foreground contamination is particularly challenging for CMB B-mode measurements. In fact, the CMB B 
mode arises from tensor perturbations (see Liddlo & Lyth 2000 for reviews), which are subdominant with respect to the 
scalar component (Spergel et al. 2003). In addition, tensor perturbations vanish on sub-degree angular scales, corresponding 
to sub-horizon scales at decoupling. On such scales, some B-mode power could be introduced by weaJi lensing (see Hu 2002 
and references therein). Anyway, the cosmological S-mode power is always expected to bo much lower than the E-mode, while 
foregrounds are expected to have approximately the same power in the two modes (Zaldarriaga 2001). 

Baccigalupi et al. (2001) estimated the power spectrum of synchrotron as derived by two main data-sets. As we already 
mentioned, on super-degree angular scales, corresponding to multipoles £ < 200, the foreground contamination is determined 
from the Brouw & Spoelstra (1976) data, covering roughly half of the sky with degree resolution. The d behavior on smaller 
angular scales has been obtained analyzing more recent data reaching a resolution of about 10' (Duncan et al. 1997, 1999, 
Uyaniker et al. 1999). These data, reaching Galactic latitudes up to 6 ~ 20°, yield a flatter slope, Ce ~ ^-(i-s-i.s) ^^^^ 
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Figure 1. Q Stokes parameter for the emission of Galactic synchrotron according to Giardino et al. (2002, left panel) and Ba<;cigalupi 
et al. (2001, right panel). The maps are in antenna temperature, at 100 GHz. 

also Tucci et al. 2000, Giardino et al. 2002) . Fosalba et al. (2002) interestingly provided evidence of a similar slope for the 
angular power spectrum of the polarization degree induced by the Galactic magnetic field as measured from starlight data. 
The synchrotron spectrum at higher frequencies was then inferred by scaling the one obtained in the radio band with a typical 
synchrotron spectral index of —2.9 (in antenna temperature). 

Giardino et al. (2002) built a full sky map of synchrotron polarized emission based on the total intensity map by Haslam et 
al. (1982), assuming a synchrotron polarized component at the theoretical maximum level of 75%, and a Gaussian distribution 
of polarization angles with a power spectrum estimated out of the high resolution radio band data (Duncan et al. 1997, 1999; 
Duncan et al. 1999). The polarization map obtained, reaching a resolution of about 10', was then scaled to higher frequencies 
by considering either a constant or a space- varying spectral index as inferred by multi-frequency radio observations. 

In this work we concentrate on low frequencies, modelling the diffuse polarized emission as a superposition of CMB and 
synchrotron. We used the synchrotron spatial template by Giardino et al. (2002), hereafter So model, as well as another 
synchrotron template, hereafter indicated as 5s, obtained by scaling the spherical harmonics coefficients of the So model to 
match the spectrum found by Baccigalupi et al. (2001). The Q Stokes parameter for the two spatial templates, in antenna 
temperature at 100 GHz, are shown in Fig. 1, plotted in a non-linear scale to highlight the behavior at high Galactic latitudes. 
Note how the contribution on smaller angular scales is larger in the So model. This is evident in Fig. 2 where we compare the 
power spectra of the So and Sb models with the CMB one, for the cosmological parameters of Eq. (1). Both models imply 
a severe contamination of the CMB E mode on large angular scales, say £ ^ 200, which remains serious even if the Galactic 
plane is cut out: cutting the region |6| < 20° decreases the Sg and Sb signals by about a factor of 10 and 3, respectively; the 
difference comes from the fact that the Sb power is concentrated more on large angular scales, which propagate well beyond 
the Galactic plane (see Fig. 2). 

In the Sg case the contamination is severe also for the first CMB acoustic oscillation in polarization, as a result of the 
enhanced power on small angular scales with respect to the Sb models (see also Fig. 1). On smaller scales both models predict 
the dominance of CMB E modes. On the other hand, CMB B modes are dominated by foreground emission, even if the region 
around the Galactic plane is cut out as we commented above. 

To obtain the synchrotron emission at different frequencies, we consider either a constant antenna temperature spectral 
index of —2.9, slightly shallower than indicated by the WMAP first year measurements (Bennett et al. 2003b), as well as a 
varying spectral index. In Fig. 3 we show the map of synchrotron spectral indices, in antenna temperature, which we adopt 
following Giardino et al. (2002). Note that this aspect is relevant especially for component separation, since all methods 
developed so far require a "rigid" frequency scaling of all the components, which means that all components should have 
separable dependencies on sky direction and frequency. Actually this requirement is hardly satisfied by real signals, and 
by synchrotron in particular. However, as we see in the next Section, FastICA results turn out to be quite stable as this 
assumption is relaxed, at least for the level of variation in Fig. 3. This make this technique very promising for application to real 
data. A more quantitative study of how the FastICA performance gets degraded when realistic signals as well instrumental 
systematics are taken into account will be carried out in a future work. 

3 COMPONENT SEPARATION FOR POLARIZED RADIATION 

Component separation has been implemented so far for the total intensity signal (see Maino et al. 2002, Stolyaxov et al. 2002, 
and references therein). In this Section we expose how we extend the ICA technique to treat polarization measurements. 
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Figure 2. E (solid) and B (dotted) angular power spectra of CMB and synchrotron polarization emission, in antenna temperature, at 
100 GHz, according to the So (left panel), and Sb (right panel) polarized synchrotron template. 




Figure 3. Map of synchrotron spectral indices (Giardino et al. 2002). 

3.1 E and B modes 

As we stressed in the previous Section, the relevant information for the CMB polarized signal can be conveniently read in 
a non-local combination of Q and U Stokes parameters, represented by the E and B modes (see Zaldarriaga & Seljak 1997 
and Kamionkowski, Kosowsky, & Stebbins 1997). There arc conceptually two ways of performing component separation in 
polarization observations. Q and U can be treated separately, i.e. performing separation for each of them independently. 
However, in the hypothesis that Q and U have the same statistical properties, separation can be conveniently performed on 
a data-set combining Q and U maps. This is surely the appropriate strategy if one is sure that the choice of polarization 
axes of the instrumental set-up does not bias the signal distribution. In general, however, it may happen that accidentally the 
instrumental polarization axes are related to the preferred directions of the underlying signal, making Q and U statistically 
different, so that merging them in a single data-set would not bo appropriate. While for the Gaussian CMB statistics we do 
not expect such occurrence, it may happen for foregrounds, expecially if separation is performed on sky patches. For example, 
the Galactic polarized signal possesses indeed large scale structures with preferred directions, such as the Galactic Plume, 
discovered by Duncan et al. (1998) in the radio band, extending up to 15 degrees across the sky and reaching high Galactic 
latitudes. Therefore, in general, the most conservative approach to component separation in polarization is to perform it for 
Q and U separately. Note that in our Galactic model no coherence in the polarization detection is present. This allow us to 
verify, in the next Section, that the results obtained by merging Q and (7 in a single data-set arc quite equivalent or more 
accurate than those obtained by treating them separately. In the following we report the relevant FastICA formalism for the 
latter case. Otherwise, when Q and U form a single template, the same formulas developed in Maino et al. (2002) do apply. 
Let these multi-frequency maps be represented by x*^ and respectively, where x is made of two indexes, labelling 
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frequencies on rows and pixel on columns. If the unknown components to be recovered from the input data scale rigidly in 
frequency, which means that each of them can be represented by a product of two functions depending on frequency and space 
separately, we can define a spatial pattern for them, which we indicate with and s^. Then we can express the inputs x*^'^ 
as 

x«'^ = A«-^Bs«'^ + n«'^ , (2) 

where the matrix A*^'^ scales the spatial patterns of the unknown components to the input frequencies, thus having a number 
of rows equal to the number of input frequencies. The instrumental noise n has same dimensions as x. The matrix B represents 
the beam smoothing operation: we recall that at the present level of architecture, an ICA based component separation requires 
to deal to maps having equal beams at all frequencies. Separation is achieved in the real space, by estimating two separation 
matrices, W*^ and W^, having a number of rows corresponding to the number of independent components and a number of 
columns equal to the frequency channels, which produce a copy of the independent components present in the input data: 

y = W X , y = W X . (3) 

All the details on the way the separation matrix for FastICA is estimated are given in Maino et al. (2002). y'^ and can 
be combined together to get the E and B modes for the independent components present in the input data (see Zaldarriaga, 
Seljak 1997 and Kamiorikowski, Kosowsky, Stebbins 1997). Note that failures in separation for even one of Q and U affect 
in general both E and B, since each of them receive contributions from both Q and U. It is possible, even in the noisy case, 
to check the quality of the resulting separation by looking at the product WA, which should be the identity in the best 
case. This means that the frequency scalings of the recovered components can be estimated: following Maino et al. (2002), 
by denoting as x'^'^ the j-th component in the data x at frequency v, it can be easily seen that the frequency scalings are 
simply the ratios of the column elements of the matrices W~^: 

^Q,u (wQ.^)-i. • ^ ^ 

However, it is important to note that even if separation goes virtually perfect, which means that WA is exactly the identity, 
Eqs. (2,3) imply that noise is transmitted to the FastICA outputs, even if it can be estimated and, to some extent, taken 
into account during the separation process. 



3.2 Instrumental noise 

Our method to deal with instrumental noise in a FASTlCA-based separation approach is described in Maino et al. (2002), 
for total intensity maps. Before starting the separation process, the noise correlation matrix, which for a Gaussian, uniformly 
distributed noise is null except for the noise variances at each frequency on the diagonal, is subtracted from the total signal 
correlation matrix; the "denoised" signal correlation matrix enters then as an input in the algorithm performing separation. 
The same is done also here, for Q and U separately. Moreover, in Maino et al. (2002) we described how to estimate the noise 
of the FastICA outputs. In a similar way, let us indicate the input noise patterns as n*^'^. Then from Eqs. (2,3) it can be 
easily seen that the noise on FastICA outputs is given by 

which moans that, if the noises on different channels are uncorrelated, and indicating as ffvj the input noise rms at frequency 

Vj, the noise rms on the i-th FastICA output is 



lY.\^r?\-?f\' ■ (6) 



Note that the above equation describes the amount of noise which is transmitted to the outputs after the separation matrix 
has been found, and not how much the separation matrix is affected by the noise. As treated in detail in Maino et al. (2002), 
if the noise correlation matrix is known, it is possible to subtract it from the signal correlation matrix, greatly reducing the 
influence of the noise on the estimation of W*^'^; however, sample variance and in general any systematics will make the 
separation matrix noisy in a way which is not accounted by Eq.(6). 

On whole sky signals, the contamination to the angular power spectrum coming from a uniformly distributed, Gaussian noise 
characterized by c is = A-Ka"^ /N, where is the pixel number on the sphere. The noise contamination to Cts on Q and 
U can therefore be estimated easily on FastICA outputs once cr^.'^ in Eq. (6) are known. Gaussianity and uniformity make 
also very easy to calculate the noise level on E and B modes, since they contribute at the same level. Thus we can estimate 
the noise contamination in the E and B channels as 

N 



^ivi - Wj/i A iTt ' \') 
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where the factor 4 is due to the normahzation according to the HEALPix scheme (version 1.10 and less), featuring conventions 
of Kamionkowski, Kosowsky, & Stebbins (1997), whereas the other common version (Zaldarriaga & Seljak 1997) would yield 
a factor 2. The quantities defined in Eq. (7) represent the average noise power, which can be simply subtracted from the 
output power spectra by virtue of the uncorrelation between noise and signal; the noise contamination is then represented by 
the power of noise fluctuations around the mean [Eq. (7)]: 



Note that the noise estimation is greatly simplified by our assumptions: a non-Gaussian and/or non-uniform noise, as well as 
a non-zero Q/U noise correlation etc. could lead for instance to non-flat noise spectra for Ces, as well as non-equal noises in E 
and B modes. However, if a good model of these effects is available, a Monte Carlo pipeline is still conceivable by calculating 
many realizations of noise to flnd the average contamination to E and B modes to be subtracted from outputs instead of the 
simple forms of Eqs. (7,8). 



4 PERFORMANCE STUDY 

In this Section we apply our approach to simulated skies to assess: (i) the ultimate capability of FastICA to clean the CMB 
maps from synchrotron in ideal noise-less conditions, (n) how the results are degraded by noise. 

4.1 Noiseless separation 

We work with angular resolution of 3'. 5, corresponding to Hgide = 1024 in an HEALPix environment (Gorski et al. 1999); 
this is enough to tost the performance of the CMB polarization reconstruction, in particular for the undamped sub-degree 
acoustic oscillations, extending up to £ ~ 2000 in Fig. 2. In all the cases we show, the computing time to achieve separation 
was of the order of a few minutes on a Pentium IV 1.8 GHz processor with 512 Mb RAM memory. We perform separation by 
considering the CMB model defined by Eq. (1) and both the So (Giardino ct al. 2002) and the Sb (Baccigalupi et al. 2001) 
model for synchrotron emission. We have considered Q and U maps both separately and combined. 

As we already mentioned, the FastICA performance turns out to be stable against relaxation of rigid frequency scalings, 
at least for the spectral index variations shown in Fig. 3. In order to illustrate quantitatively this point, we compare the quality 
of the CMB reconstruction assuming either constant and varying synchrotron spectral indices, (3. In Fig. 4 we plot the original 
(dotted) and reconstructed (solid) ' s for CMB, in the case of the Sb foreground model with constant (left-hand panels) or 
spatially varying /3. The upper (lower) panels refer to the two frequency combinations 70, 100 (30, 44) GHz channels. Figure 5 
shows the results of the same analisys, but using the So model. 

As it can be seen, the CMB signal is well reconstructed on all relevant scales, down to the pixel size. The same is true for 
the synchrotron emission, not shown. Percent precision in frequency scaling recovery for CMB and synchrotron is achieved (sec 
Table 1). As we stressed in the previous Section, the precision on frequency scaling recovery corresponds to the precision on 
the estimation of the elements in the inverse of the matrices W*^ and W^. Remarkably, FastICA is able to recover the CMB 
B modes on all the relevant angular scales, even if they are largely subdominant with respect to the foreground emission, as 
it can be seen in Fig. 2. This is due to two main reasons: the difference of the underlying statistics describing the distribution 
of CMB and foreground emission, and the high angular resolution of templates (3'. 5). Such resolution allows the algorithm to 
converge close to the right solution by exploiting the wealth of statistical information contained in the maps. Note also that 
in the noiseless case with constant synchrotron spectral index, the CMB power spectrum is reconstructed at the same good 
level both for the 100, 70 GHz and for 30, 44 GHz channel combinations, although in this frequency range the synchrotron 
emission changes amplitude by a factor of about 10. Indeed, by comparing top and bottom left panels of Figs. 4 and 5, one 
can note that there is only a minimum difference between the B spectra at 44 GHz and at 100 GHz, arising at high i, while 
the E spectra exhibit no appreciable difference at all. 

In the case of spatially varying spectral index (right-hand panels of Figs. 4 and 5) a rigorous component separation is 
virtually impossible, since the basic assumptions of rigid frequency scaling is badly violated. However FastICA is able to 
approach convergence by estimating a sort of "mean" foreground emission, scaling roughly with the mean value of the spectral 
index distribution. Some residual synchrotron contamination of the CMB reconstructed maps cannot be avoided, however. 
This residual is proportional to the difference between the "true" synchrotron emission and that corresponding to the "mean" 
spectral index and is thus less relevant at the higher frequencies, where synchrotron emission is weaker. As shown by the upper 
right-hand panels of Figs. 4 and 5, when the 70-100 GHz combination is used, the power spectrum of the CMB i?-mode is 
still well reconstructed on all scales, and even that of the CMB B-mode is recovered at least up to £ ~ 100. On smaller scales 
synchrotron contamination of the B-mode is strong in the So, but not in the Sb (at least up to ^ ~ 1000) case. As expected, 
the separation quality degrades substantially if the 30-44 GHz combination is used (right-hand, bottom panels of Figs. 4 and 




(8) 
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Sg with constant spectral index 




Sg witti varying spectral index 




Sg wltti constant spectral index 




Sg witti varying spectral Index 




Figure 4. Original (dotted) and reconstructed (solid) power spectra for E and B CMB modes in the case of the Sb model, in the 
absence of noise. Left panels: constant spectral index. Right panels: space varying spectral index. Upper panels are from inputs at 70 and 
100 GHz channels, while bottom panels corresponds to inputs at 30 and 44 GHz. The outputs are conventionally plotted at the highest 
input frequency. 



5; see also Table 1 where the quoted error on frequency scaling for varying spectral index is the percentage difference between 
the average values < (i'i/i'2Y > computed on input and reconstructed synchrotron maps). 

To compare the FastICA performances when Q and U maps are dealt with separately or together (case QU), we have 
carried out a Monte Carlo chain on the CMB realizations, referring to the 70 and 100 GHz channels, and computed the 
rms error on the CMB frequency scaling reconstruction, aq, au, and aqu- The results are shown in Table 2. Again the 
reconstruction is better when the weaker synchrotron model Sb is considered. The slight difference between aq and au is 
probably due to the particular realization of the synchrotron model we have used (not changed through this Monte Carlo 
chain). The fact that the difference is present for both the Sb and So models is not surprising because the two models have 
different power as a function of angular scales but have the same distribution of polarization angles. The reason why we didn't 
vary the foregrounds templates in our chain is the present poor knowledge of the underlying signal statistics. 

The separation precision when Q and U axe considered together is equivalent or better than if they are treated separately, 
as expected since the statistical information in the maps which are processed by FastICA is greater. On the other hand, the 
fact that aqu is so close to aq and au clearly indicates that for a pixel size of about 3'. 5 the statistical information in the 
maps is such that the results are not greatly improved if the pixel number is doubled. In the following we just consider the 
most general case in which Q and U maps are considered separately. 
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Figure 5. Original (dotted) and reconstructed (solid) power spectra for E and B CMB modes in the case of Sq model, in the absence 
of noise. Left panels: constant spectral index. Right panels: space varying spectral index. Upper panels are from inputs at 70 and 100 
GHz channels, while bottom panels corresponds to inputs at 30 and 44 GHz. The outputs are conventionally plotted at the highest input 
frequency 



4.2 The effect of noise 

To study the effect of the noise on FastICA component separation we use a map resolution of about 7', corresponding to 
nside = 512 in an HEALPix scheme (Gorski et al. 1999). At this resolution the all sky separation runs take a few seconds. 
Moreover we consider only the combination of 70 and 100 GHz channels, and a space varying synchrotron spectral index. 
We give the results for one particular noise realization and then we show that the quoted results are representative of the 
typical FastICA performance within the present assumptions. Moreover, we investigate how the foreground emission affects 
the recovered CMB map. The noise is assumed to be Gaussian and uniformly distributed, with rms parameterized with the 
signal to noise ratio, S/N, where S stays for CMB. As we already stressed, noise is either subtracted during the separation 
process, and on the reconstructed Ce, according to the estimate in Eq, (7). The results are affected by the residual noise 
fluctuations, with power given by Eq. (8). 

As expected, the noise primarily affects the reconstruction of CMB B mode. In Fig. 6 we plot the reconstructed and 
original CMB E- and S-mode power spectra, for So (left) and Sb (right) foreground emission, in the case S/N = 2. With 
this level of noise, separation is still successful: the _B-mode power spectrum comes out very well, while that of B-mode is 
well reconstructed up to the characteristic peak at ^ ~ 100. Table 3 shows that the error on the frequency scaling recovery, 
for CMB, remains, in the noisy case, at the percent level both for Q and U. 

If the S/N ratio is decreased, B modes get quickly lost, while still the algorithm is successful in recovering the E power 
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Table 1. Percentage errors on frequency scalings reconstruction in the noiseless case. 



Synchrotron S(y model 


const. /3 




variable f3 


Frequencies used (GHz) 


30, 44 


70, 100 




30, 44 


70, 100 


Q synch. 


1.48 xlO-2 


0.11 




21.73 


7.55 


U synch. 


5.71 xlO-3 


4.56 xlO" 


.2 


21.73 


7.55 




0.31 


2.69 xlO" 


-2 


207.36 


1.18 


U CMB 


0.69 


5.83 xlO" 


-2 


222.12 


1.21 


Synchrotron Sb model 


const. /3 




variable f3 


Frequencies used (GHz) 


:?(). 11 


70, 100 




:',(). 11 


70, 100 


Q synch. 


0.18 


1.46 




19.19 


1.04 


U synch. 


0.137 


1.10 




19.73 


1.04 


Q CMB 


1.01 


8.79 xlO" 


-2 


6.15 


9.62 xlO"^ 


U CMB 


0.1448 


1.24 xlO" 


-2 


9.028 


0.107 



Table 2. Percentage rms for CMB frequency scaling reconstruction resulting from a Monte Carlo chain of FastICA applied to 50 
different CMB realizations at 70 and 100 GHz; aq and ctj/ are obtained by treating Q and U separately, while uqu is the result when 
they are considered as a single array; the la error on the parameter estimation, assuming Gaussianity, is also indicated. 









o-(7 




Sy nchr ot ron 


Sg model 


1.23 ±0.17 


1.17 ±0.17 


1.21 ±0.17 


Sy nchr ot ron 


Sb model 


0.104 ±0.015 


0.113 ±0.016 


0.102 ±0.014 



spectrum for S/N ^ 0.2 (see Fig. 7). The algorithm starts faiUng at low multipoles, say £ ^ 100, where synchrotron dominates 
over the CMB both in the So and Sb cases. The results in Fig. 7 for the Sg case are averages over 8 multipoles, to avoid 
excessive oscillations of the recovered spectrum. For the S/N ratios in this figure, the B-modes are lost on all scales. Since 
the Sg model has a higher amplitude, FastICA is able to catch up the statistics more efficiently than for Sb, thus being able 
to work with a lower S/N . Table 3 shows the degradation of the separation matrix for the S/N ratios of Fig. 7, compared to 
the case with S/N = 2. 

We stress that the noise levels quoted here are not the maximum which the algorithm can support. The quality of the 
separation depends on the noise level as well as on the number of channels considered; adding more channels, while keeping 
constant the number of components to recover, generally improves the statistical sample with which FastICA deals and so 
the quality of the reconstruction as well as the amount of noise supported. In the next Section we show an example where a 
satisfactory separation can be obtained with higher noise by considering a combination of three frequency channels. 

Table 3. Percentage errors on CMB frequency scalings reconstruction in the noisy case by considering the 70 and 100 GHz channels. 



Synchrotron 


Sg model 


S/N = 2 


S/N = 0.2 


Q CMB 




1.04 


13 


U CAIL! 




O.OG 


0.S9 


Sy nchr ot ron 


Sb model 


S/N = 2 


S/N = 0.5 


Q CMB 




0.31 


4.26 


U CMB 




0.34 


0.35 
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Sq with varying spectral index Sg witti varying spectral index 





Figure 6. Original (dotted) and reconstructed (solid) 100 GHz power spectra for E (top) and B (bottom) CMB modes for the Sc (left) 
and Sb (right) cases, assuming S/N=2 and considering the 70 and 100 GHz channels; the synchrotron spectral index is space varying. 



We now investigate to what extent the results quoted here are representative of the typical FastICA performance and 
study how the foreground emission biases the CMB maps recovered by FastICA. To this end we have performed a Monte 
Carlo chain of separation runs, building for each of them a map of residuals by subtracting the input CMB template from the 
recovered one and studying the ensemble of those residual maps. 

The residuals in the noiseless case are just a copy of the foreground emission. Their amplitude is greatly reduced with 
respect to the true foreground amplitude, in proportion to the accuracy of the recovered separation matrix. Since the latter 
accuracy is at the level of percent of better, as it can be seen in Table 1, the residual foreground emission in the CMB recovered 
map is roughly the true one divided by 100. In terms of the angular power spectrum, see e.g. Fig. 2, the residual foreground 
contamination to the CMB recovered power spectrum is roughly a factor 10* less than the true one. 

In the of noisy separation a key feature is that at the present level of architecture, the FastICA outputs are just a linear 
combination of the input channels. Thus, even if the separation goes perfect, the noise is present in the output just as the 
same linear combination of the input noise templates. Note however that this docs not mean that the noise is transmitted 
linearly to the outputs. The way the separation matrix is found depends non-linearly on the input data including the noise. In 
other words, the noise affects directly the estimation of the separation matrix, as we explained in Section 3.2. Eq. (6) describes 
only the amount of noise which affects the outputs after the separation matrix is found. As we shall see in a moment, at least 
in the case S/N = 2 the main effect of the noise is the one given by Eq. (6), dominant over the error induced by the noise on 
the separation matrix estimation. Moreover, the noise in the outputs reflects the input noise statistics, which is Gaussian and 
uniformly dsitributed in the sky. As we shall see now, this is verified if the foreground contamination is the stronger So one. 
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Table 4. Statistics of CMB residuals in Kelvin and percentage errors on frequency scaling reconstruction for the ease Sg, S/N = 2, 70 
and 100 GHz channels, on 50 different noise and CMB realizations. The la error on the parameter estimation, assuming Gaussianity, is 
also indicated. 



Stokes parameter < r > < >q a 

Q (3.94 ± 0.56) • 10"® (2.20 ± 0.31) • 10-^ 1.14 ±0.16 

(/ ( 1.03 ± 0.57) • 10"^ (2.20 ± 0..31) • 10"" 0,71 ± 0.10 



The results presented in Table 4 show the ensemble average of the mean of the residuals < f > together with its Gaussian 
expectation < >g^^ and the mean rms error on the CMB frequency scaling recovery, a. The most important feature is 
that a non-zero mean value, at almost 10a with respect to its Gaussian expectation, is detected. This is the only foreground 
contamination we find in the residuals. Note that the separation matrix precision recovery is at the percent level. That 
means that the present amount of noise does not affect significantly the accuracy of the separation process. Of course, if the 
noise is increased, the separation matrix estimation starts to be affected and eventually the foreground residual in the CMB 
reconstructed map will be relevant. 

We made a further check by verifying that the residuals obey a Gaussian statistics with rms given by Eq. 6 on all 
Galactic latitudes. We constructed a map having in each pixel the variance built out of the 50 residual maps in our Monte 
Carlo chain. In Fig. 8 we show the rms of such map, plus/minus the standard deviation, calculated on rings with constant 
latitude with width equal to 1 degree. Together with the curves built out of our Monte Carlo chain, we report the theoretical 
values according to a Gaussian statistics, i.e. the average given by Eq. (6) equal to (3.90 x 10~^K)^, and the standard deviation 
over N = 50 samples, given by y/2{N — 1)/N^ ■ (3.90 x 10~^K)^. The agreement demonstrates that the Gaussian expectation 
is satisfied at all latitudes, expecially at the lowest, where the foreground contamination is expected to be maximum. Note 
that the fluctuations around the Gaussian theoretical levels are larger near the poles because of the enhanced sample variance. 

We conclude that, within the present assumptions, for a successful separation the residual foreground contamination in 
the recovered CMB map is subdominant with respect to the noise. On the other hand, further tests are needed to check this 
result against a more realistic noise model, featuring the most important systematic effects like a non-uniform sky distribution, 
the presence of non-Gaussian features etc. 



5 AN APPLICATION TO PLANCK 

In this Section, we study how FastICA behaves in conditions corresponding to the instrumental capabilities of Planck. 
While this work was being completed, the polarization capabilities at 100 GHz were lost because of a funding problem of 
the Low Frequency Instrument (LFI), but that capability could be restored if the 100 GHz channel of the High Frequency 
Instrument (HFI) is upgraded, as it is presently under discussion. The Planck polarization sensitivity in all its channels has 
a crucial importance, and it is our intention to support this issue. Thus we work assuming Planck polarization sensitivity at 
100 GHz, highlighting the fact that of our results have been obtained under this assumption. At 30, 44, 70 and 100 GHz, the 
Planck beams have full width half maximum (FWHM) of 33, 23, 14 and 10 arcmin, respectively. We study the FastICA 
effectiveness in recovering E, B and TE modes, separately. We adopted, for polarization, the nominal noise level for total 
intensity measurements increased by a factor (note also that due to the 1.10 and lower HEALPix version convention to 
normalize Q and U following the prescription by Kamionkowski, Kosowsky and Stebbins 1997, a further \/2 has to be taken 
into account when generating Q and U maps out of a given power in E and B). We neglect all instrumental systematics in 
this work. The Planck instrumental features assumed here, with noise rms in antenna temperature calculated for a pixel 
size of about 3.52 arcmin corresponding to nside — 1024 in the HEALPix scheme, arc summarized in Table 5. By looking at 
the numbers, it can be immediately realized that the level of noise is sensibly higher than the one considered in the previous 
Section, so that the same method would not work in this case and an improved analysis, involving more channels as described 
below, is necessary. 

5.1 E mode 

Due to the high noise level, we found convenient to include in the analysis the lower frequency channels together with those 
at 70 and 100 GHz. Since the FastICA algorithm is unable to deal with channels having different FWHM, as in Maino et 
al. (2002), we had to degrade the maps, containing both signal and noise, to the worst resolution in the channels considered. 
However, a satisfactory recovery of the CMB E modes, extending on all scales up to the instrument best resolution, is still 
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Figure 7. Input (dotted) and reconstructed (solid) 100 GHz power spectra for E (top) and B (bottom) CMB modes for tlie Sq (left, 
assuming S/N = 0.2, averaged over 8 multipole intervals) and Sb (right, assuming S/N = 0.5) cases, using the 70 and 100 GHz channels. 
The synchrotron spectral index is space varying. 



Table 5. PLANCK-polarization performance assumed in this work. 



Frequency (GHz, LFI) 


30 


44 


70 


100 


FWHM (arcmin.) 


33 


23 


14 


10 


noise rms for 3.52' pixels (/uK) 


80 


81 


69 


52 



possible by making use of the different angular scale properties of both synchrotron and CMB. Indeed, as it can be seen in 
Fig. 2, the Galaxy is likely to be a substantial contaminant on low multipoles, say £ ^ 200. 

For the present application, we found convenient to use a combination of three Planck channels, 44, 70, 100 GHz for 
the Sg and 30, 70, 100 GHz for the Sb models, respectively. The reason of the difference is that the Sb contamination is 
weaker, and the 30 GHz channel is necessary for FastICA to catch synchrotron with enough accuracy. Including a fourth 
channel does not imply relevant improvements. The maps, including signals properly smoothed and noise according to the 
Table 5, were simulated at 3.52' resolution, corresponding to Uside = 1024. Higher frequency maps were then smoothed to the 
FWHM of the lowest frequency channel and then re-gridded to Uside = 128, corresponding to a pixel size of about 28' and to 
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Galactic latitude [degrees] 

Figure 8. Latitude analysis of the variance map calculated out of the residual maps in our Monte Carlo chain, in the case Sq, S/N = 2, 
70 and 100 GHz channels. The solid line is the average at the corresponding latitude, while the dashed curve represent the average 
plus/minus the standard deviation. The dotted lines are derived assuming Gaussian statistics. 





Figure 9. Original (dotted) and reconstructed (solid) CMB obtained by applying the FastICA algorithm to the combination of 44, 
70, and 100 GHz channels for the Sq synchrotron case (left) and of 30, 70, and 100 GHz channels for the Sb model (right). 



a maximum multipolc £ ~ 400. In all the cases shown, the spectral index for synchrotron has been considered variable. Fig. 9 
shows the resulting CMB E mode power spectrum after separation, for the So (left) and Sb (right) synchrotron models. An 
average every 4 (left) and 3 (right) coefficients was applied to eliminate fluctuations going negative on the lowest signal part 
at / ~ 10. The agreement between the original spectrum and the reconstructed one is good on all the scales probed at the 
present resolution, up to £ ~ 400. The re-ionization bump is clearly visible as well as the first polarization acoustic oscillation 
at £ ~ 100. Moreover, there is no evident difference in the quality of the reconstruction between the two synchrotron models 
adopted. 

Let us turn now to the degree and sub-degree angular scales, £ ^ 200. As we already stressed, the Galaxy is expected to 
yield approximately equal power on E and B modes (see Fig. 2). On the other hand, CMB E and B modes are dramatically 
different on sub-degree angular scales. Summarizing, on. £ ^ 200, we expect to have 



^ECMB „ ^BGMB 



(9) 



Therefore, on £ ^ 200 where the CMB contamination from synchrotron is expected to be irrelevant, the power spectrum of 
the CMB E modes can be estimated by simply subtracting, together with the noise, the B power as 



CMB+Gal 



-c, 



B CMB+Gal 



(10) 



where the total map CMB+Gal is used without any separation procedure. In other words, there is no need to perform 
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Figure 10. Input (dotted) and reconstructed (solid) CMB obtained by subtracting tlie expected level of noise as well as the 
synchrotron contaminations Sq (left) and Sb (right) assumed to be matched by the B-mode map. The adopted instrumental capabilities 
are those of the Planck-HFI channel at 100 GHz, assumed to be polajrization-sensitive. 

separation to get the CMB E modes at high multipoles, because they are simply obtained by subtracting the B modes of the 
sky maps, since the latter arc dominated by synchrotron which has almost equal power on E and B modes. Fig. 10 shows 
the results of this technique applied to the Planck-HFI 100 GHz channel, assumed to have polarization capabilities, for both 
the So and Sb synchrotron models. Residual fluctuations are higher in the So case since the synchrotron contamination 
is stronger. In both cases, CMB E modes are successfully recovered in the whole interval 100 ^ I ^ 1000. It has also to be 
noted that the same subtraction technique would not help on the lower multipoles considered before, since in that case the 
foreground contamination is so strong that the tiny fluctuations making Galactic E and B modes different are likely to hide 
the CMB signal anyway. 

Our results here can be summarized as follows: in the case of Planck capabilities, the FastICA technique makes 
it possible to remove substantially the foreground contamination in the regions in which that is expected to be relevant. 
Planck is likely to measure the CMB E modes over all multipoles up to £ ~ 1000. 

5.2 B mode 

In Fig. 11 the B-mode power spectrum after FastICA separation described in the previous section is shown for the So (left) 
and Sb (right) synchrotron models. An average over 13 multipoles has been applied to both cases in order to avoid fluctuations 
going negative. The reconstructed signal approaches the original one at very low multipoles, say £ ^ 5. At higher multipoles, 
where the B signal is generated by gravitational waves, the overall amplitude appears to be recovered, even if with major 
contaminations especially in the region where the signal is low, i.e. right between the re-ionization bump and the raise toward 
the peak at £ ~ 100. Needless to say, such contaminations are due to a residual foreground emission. 

In the insets of Fig. 11 we show (data points) the recovered B-mode power spectrum in the range between 30 < i < 120, 
averaged over 20 multipoles, with error bars given by Eq. (8) . Even if the contamination is substantial, especially for £ > 100 
and for the Sg case, the results show a sign of the characteristic rise of the spectrum due to cosmological gravitational waves. 

Concluding, our results indicate that the FastICA technique is able to remove substantially the foreground contamination 
of the B mode, up to the peak at £ ^ 100 if the tensor to scalar perturbation ratio is at least 30%. 

5.3 TE mode 

While the cosmological TE power spectrum is substantially stronger than that of any other polarized CMB mode, the opposite 
should happen in the case of foregrounds. On degree and sub-dcgrcc angular scales, a measure of the synchrotron TE power 
spectrum can be achieved in the radio band. In the Parkes data at 1.4 GHz, Uyaniker et al. (1999) were able to isolate a 
region exhibiting pretty low rotation measures, named "fan region" , which is therefore expected to be only weakly affected 
by Faraday depolarization. This and other regions from the existing surveys in the radio band were used to predict the 
synchrotron power for the Sb scenario (Baccigalupi et al. 2001). 

In Fig. 12 we show the T, E, B, and TE power spectra for the fan region. Total intensity anisotropies axe represented 
by the upper curve (solid). E and B modes (light lines) have very similar behavior. The TE mode (heavy solid line) is the 
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Figure 11. Original (dotted) and reconstructed (solid) CMB obtained by applying the FastICA algorithm to the combination of 
Planck 44, 70, and 100 GHz channels for the Sq synchrotron case (left) and to the combination of 30, 70, and 100 GHz channels for 
the Sb model (right). The data points in the insets show the recovered B-mode power spectrum in the range 30 < £ < 120 averaged 
over 20 multipoles; the error bars are given by Eq. (8). 




Figure 12. Power spectra of synchrotron T (solid linc),i? and B (light lines), and TE (heavy solid line) modes in the fan region at 
medium Galactic latitudes, at 1.4 GHz (Uyaniker et al. 1999). 



weakest and, as it can be easily seen by scaling the TE amplitude in Fig. 12 with the typical spectra index for synchrotron, 

it is markedly below the expected cosmological TE signal at CMB frequencies. Both synchrotron models So and Sb, are 
consistent with this result as illustrated by in Fig. 13. It is straightforward to check that our models for the synchrotron 
emission have a TE power spectrum not fax from the one in Fig. 12, when scaled to the appropriate frequency. 

From the point of view of CMB observations, this moans that, if the synchrotron contamination at microwave frequencies 
is well represented by its signal in the radio band, at least on degree and sub-degree angular scales the contamination from 
synchrotron is almost absent due to the change in the magnetic field orientation along the line of sight. On the other hand, on 
larger scales, as it can be seen in Fig. 13, the contamination could be relevant both in the So and Sb cases and we perform 
component separation as described in Section 5.1 for the Planck case. In Fig. 14 we show the recovery of the CMB TE 
mode, obtained by combining the templates of Q and U maps obtained after FastICA application as in Section 5.1, with 
the CMB T template obtained, still with FastICA based component separation strategy, in Maino et al. (2002). Oscillations 
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Figure 13. Cj^^ of CMB (dotted) compajred with that of synchrotron according to the So (solid, left) and Sb (solid, right) models, at 
100 GHz. 




Figure 14. Cj'^ of original (dotted) and recovered (solid) CMB emission obtained with FastICA applied to simulated Planck maps, 
by considering the Sq (left) and the Sb (right) foreground model, respectively. Results are shown at 100 GHz. 

due to residual noise are visible in the recovered Cf^. However, as in the case of the E-mode, the procedure was successful 
in substantially removing the contamination. 

Both in the Sg and in the Sb case the synchrotron contamination is almost absent in the acoustic oscillation region of the 
spectrum, as it is evident again from Fig. 13; neglecting it wc get the results shown in Fig. 15. The combination of Planck 
angular resolution and sensitivity allows the recovery of the TE power spectrum up to £ ~ 1200, corresponding roughly to 
the seventh CMB acoustic oscillation. 



6 CONCLUDING REMARKS 

Forthcoming experiments are expected to measure CMB polarization^. The first detections have been obtained on pure 
polarization (Kovac et al. 2002), as well on its correlation with total intensity CMB anisotropics, by the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite'^ . 

The foreground contamination is mildly known for total intensity measurements, and poorly known for polarization (see 

® see lambda.gsfc.nasa.gov/ for a collection of presently operating and future CMB experiments 
^ map.gsfc.nasa.gov/ 
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Figure 15. Cj^ of original (dotted) and recovered (solid) CMB emission adopting the Sq (left) and the Sb (right) synchrotron model, 
considering the Planck performances at 100 GHz. 

De Zotti 2002 and references therein). It is therefore crucial to develop data analysis tools able to clean the polarized CMB 
signal from foreground emission by exploiting the minimum number of a priori assumptions. In this work, we implemented the 
Fast Independent Component Analysis technique in an astrophysical context (FastIC A, see Amari, Chichocki 1998, Hyvarinen 
1999, Baccigalupi et al. 2000, Maino et al. 2002) for blind component separation to deal with astrophysical polarized radiation. 

In our scheme, component separation is performed both on the Stokes parameters Q and U maps independently and by 
joining them in a single dataset. E and B modes, coding CMB physical content in the most suitable way (see Zaldarriaga, 
Seljak 1997, Kamionkowski, Kosowsky, Stebbins 1997), are then built out of the separation outputs. We described how to 
estimate the noise on FastICA outputs, on Q and U as well as on iJ and B. 

We tested this strategy on simulated polarization microwave all sky maps containing a mixture of CMB and Galactic 
synchrotron. CMB is modelled close to the current best fit (Spergel et al. 2003), with a component of cosmological gravitational 
waves at the 30% level with respect to density perturbations. We also included re-ionization, although with an optical depth 
lower than indicated by the WMAP results (Bennett et al. 2003a) since they came while this work Weis being completed, but 
consistent with the Gunn-Peterson measurements by Becker et al. (2001). Galactic synchrotron was modelled with the two 
existing templates by Giardino et al. (2002) and Baccigalupi et al. (2001). These models yield approximately equal power on 
angular scales above the degree, dominating over the expected CMB power. On sub-degree angular scales, the Giardino et 
al. (2002) model predicts an higher power, but still subdominant compared to the CMB E mode acoustic oscillations. Note 
that at microwave frequencies, the fluctuations at high multipoles {£ ^ 1000), corresponding to a few arcmin angular scales, 
are likely dominated by compact or flat spectrum radio sources (Baccigalupi et al. 2002b, Mesa et al. 2002). Their signal is 
included in the maps used to estimate the synchrotron power spectrum. 

We studied in detail the limiting performance in the noiseless case, as well as the degradation induced by a Gaussian, 
uniformly distributed noise, by considering two frequency combinations: 30 + 44 GHz and 70 + 100 GHz. In the noiseless 
case, the algorithm is able to recover CMB E and B modes on all the relevant scales. In particular, this result is stable against 
the space variations of the synchrotron spectral index indicated by the existing data. In this case, FastICA is able to converge 
to an average synchrotron component, characterized by a "mean" spectral index across the sky, and to remove it efficiently 
from the map. The output CMB map, also containing residual synchrotron due to its space varying spectral index, is mostly 
good as far as the frequencies considered are those where the synchrotron contamination is weaker. 

By switching on the noise we found that separation, at least for what concerns the CMB E mode, is still satisfactory for 
noise exceeding the CMB but not the foreground emission. The reason is that in these conditions the algorithm is still able 
to catch and remove the synchrotron component efficiently. We implemented a Monte Carlo chain varying the CMB and the 
noise realizations in order to show that the performance quoted above is typical and does not depend on the particular case 
studied. Moreover, we studied how the foreground emission biases the recovered CMB map, by computing maps of residuals, 
i.e. subtracting the true CMB map out of the recovered one. In the noiseless case, the residual is just a copy of the foreground 
emission, with amplitude decreased proportionally to the accuracy of the separation matrix. In the noisy case, for interesting 
noise amplitudes the residual maps are dominated by the noise in the input data, linearly mixed with the separation matrix. 
The situation is obviously worse for the weaker CMB B mode. 

We applied these tools making reference to the Planck polarization capabilities, in terms of frequencies, angular resolution 
and noise, to provide a first example of how the FastICA technique could be relevant for high precision large polarization 
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data-sets. We addressed separately the analysis of the CMB E, B, and TE modes. While this work was being completed, 
the Low Frequency Instrument (LFI) lost its 100 GHz channel, having polarization sensitivity. However polarimetry at this 
frequency could bo restored if the 100 GHz channel of the High Frequency Instrument (HFI) is upgraded, as is presently under 
discussion. Due to the scientific content of the CMB polarization signal, the Planck polarization sensitivity deserves a great 
attention. Within our context here, it is our intention to support the importance of having polarization capabilities in all the 
cosmological channels of Planck, and in particular at 100 GHz. Our results have been obtained under this assumption. 

To improve the signal statistics, we found convenient to consider at least three frequency channels in the separation 
procedure, including the ones where the CMB is strongest, 70 and 100 GHz, plus one out of the two lower frequency channels, 
at 30 and 44 GHz. Since the latter have lower resolution wo had to degrade the higher frequency maps since the present 
FastICA architecture cannot deal with maps having different resolutions. CMB E and TE modes were accurately recovered 
for both the synchrotron models considered. The S-mode power spectrum is recovered on very large angular scales in the 
presence of a conspicuous re-ionization bump. On smaller scales, where the B-mode power mainly comes from cosmological 
gravitational waves, the recovery is only marginal for a 30% tensor to scalar perturbation ratio. 

On the sub-degree angular scales the contamination from synchrotron is almost irrelevant according to both models 
(Giardino ct al. 2002, Baccigalupi et al. 2001). Moreover, it is expected that Galactic E and B modes have approximately 
the same power (Zaldarriaga 2001), while for CMB the latter are severely damped down since they are associated with vector 
and tensor perturbations, vanishing on sub-horizon scales at decoupling corresponding to a degree or less in the sky (see Hu 
et al. 1999). This argument holds also if the B mode power is enhanced by weak lensing effects from matter structures along 
the line of sight (see Hu 2002 and references therein). Therefore, on these scales, we expect the E power spectrum to be a 
sum of Galactic and CMB contributions, while the B power comes essentially from foregrounds only. In these conditions, the 
CMB E power spectrum is recovered by simply subtracting the B power spectrum. 

We also estimated the TE contamination from synchrotron to be irrelevant for CMB, because of the strength of the CMB 
TE component due to the intrinsic correlation between scalar and quadrupole modes exciting E polarization. By applying 
these considerations on sub-degree angular scales, as well as the results of the FastICA procedure described above on larger 
scales, we show how the Planck instrument is capable of recovering the CMB E and TE spectra on all scales down to 
the instrumental resolution, corresponding to a few arc-minutes scales. In terms of multipoles, the E and TE angular power 
spectra are recovered up to 1000 and 1200, respectively. 

Summarizing, we found that the FastICA algorithm, when applied to a PLANCK-like experiment, could be able to 
substantially clean the foreground contamination on the relevant multipoles, corresponding to degree angular scales and 
above. Since the foreground contamination on sub-degree angular scales is expected to be subdominant, the CMB TE and E 
modes are recovered on all scales extending from the whole sky to a few arc-minutes. In particular, the FastICA algorithm can 
clean the B-mode power spectrum up to the peak due to primordial gravitational waves if the cosmological tensor amplitude 
is at least 30% of the scalar one. In particular, we find that on large angular scales, of a degree and more, foreground 
contamination is expected to be severe and the known blind component separation techniques are able to efHciently clean the 
map from such contamination, as it is presently known or predicted. 

Still, despite these good results, the main limitation of the present approach is the neglect of any instrumental systematics. 
While it is important to assess the performance of a given data analysis tool in the presence of the nominal instrumental 
features, as we do here, a crucial test is checking the stability of such tool with respect to relaxation of the assumptions 
regarding the most common sources of systematic errors, like beam asymmetry, non-uniform and/or non-Gaussian noise 
distribution etc., as well as the idealized behavior of the signals to recover. In this work we had a good hint about the second 
aspect, since we showed that FastICA is stable against relaxation of the assumption, common to all component separation 
algorithms developed so far, abotit the separability between space and frequency dependence for all the signal to recover. In 
a forthcoming work we will investigate how ICA based algorithms for blind component separation deal with maps affected by 
the most important systematics errors. 
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